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ABSTRACT 

\ We report the discovery of a large-scale structure of Lyman a emitters (LAEs) at z = 4.86 based on wide- 

st ■ field imaging with the prime-focus camera (Suprime-Cam) on the Subaru telescope. We observed a 25' x 45' 

^ \ area of the Subaru Deep Field in a narrow band (NB711, Ac = 7126A and FWHM= 73A) together with R and 

/'. We isolate from these data 43 LAE candidates down to NB711 = 25.5 mag using color criteria. Follow-up 
' spectroscopy of five candidates suggests the contamination by low-z objects to be ^ 20%. We find that the LAE 

candidates are clustered in an elongated region on the sky of 20 Mpc in width and 50 Mpc in length at z = 4.86, 
which is comparable in size to present-day large-scale structures (we adopt Hq = 70 km s"' Mpc~', VLq = 0.3, and 
' Ao = 0.7). This elongated region includes a circular region of 12 Mpc radius of higher surface overdensity {Sy. = 2), 

which may be the progenitor of a cluster of galaxies. Assuming this circular region to be a sphere with a spatial 
overdensity of 2, we compare our observation with predictions by Cold Dark Matter models. We find that an 
" sj \ ^0 = 0.3 flat model with erg = 0.9 predicts the number of such spheres consistent with the observed number (one 

. sphere in our survey volume) if the bias parameter of LAEs is ~ 6. This value suggests that the typical mass 

' of dark haloes hosting LAEs at z — 5 is of the order of IO'^Mq. Such a large mass poses an interesting question 

\ about the nature of LAEs. 

CSJ ' Subject headings: cosmology: observations — cosmology: early universe — cosmology: large-scale structure of 

\ universe — galaxies: high-redshift — galaxies: evolution — galaxies: photometry 

cn . 
o ■ 

. L INTRODUCTION cluster of M 10'^Mq with a size of 2.7 x 1.8 Mpc around a 

' rru c 4.- J w c ii. ^ lit- f radio galaxy at z = 4. 1 through a narrow-band survcy. 

O . The formation and evolution of the large-scale clustering of °. 

Jr^, , . ... ^ , . ^ , . To place stringent constraints on models of structure for- 

I ' galaxies, as seen in the present-day universe, are central issues . ^ . . . , , , , . 

O- • , f 1 1 » • ■ A- ^ ^ mation, it is crucial to observe galaxy clustering over a wide 

. in cosmology. Observations of galaxy clustering in the distant V , . , t . • r 

?H • . . , . .u- • • 1 J- u 1 range of scales at various epochs. In this Leffer, we report on 

> ■ universe give us great clues to this issue, including how galax- ° , , . Ji, ^ . „ . 

fj^i , ■ % j-.u ji ji T-j. a narrow-band survey in a 25 x 45 field for Lyman a Emitters 

zA les were formed in the underlying dark matter To date, many ...t-s a r,^ ^ ■ . ..w 

Cu „„^, , 11.- c ^ (LAEs) at z = 4.86. Our main goal IS to search at high redshifts 

• . efforts have been made to detect large-scale clustering of galax- ). , ■ r . , , . 

.J,, n J, , 1 ° for the progenitors of present-day large-scale structures like the 

les at high redshifts. Probable detections of galaxy clusters on ^ „ n . • 

® 6 J' Qjgjjj yf^^ii (Geller & Huchra 1989), in order to trace the evo- 



>< 



Megaparsec scales have been reported (e.g., Giavalisco et al. ,. ^ ,,- 
infi^ n 11 » 1 inn/; t f i i nn<; iv^ 11 lution of Structures on scales of tens of Mpc. We also aim at 

1994; Pascarefle et al. 1996a; Le Fevre et al. 1996; Malkan, , . , , , ^ .r , ■ . 



OBSERVATIONS 



S ' Teplilz, & McLean 1996; Francis et al. 1996; Keel et'al. 1999! "^^'''^^ clusters at z ^ 5. Unless otherwise noted, we 

Campos et al. 1999; Pentericci et al. 2000). The most re- ^^'^70! km s > M c ^ 
markable observation is the discovery of a high concentration " ^" P 

of Lyman-break Galaxies (LBGs) at z ~ 3 with a size of at least 
11' X 8' (or 21 X 15 Mpc) by Steidel et al. (1998), who ar- 
gue that LBGs must be very biased tracers of mass if such a 2.1. Imaging 

structure is consistent with Cold Dark Matter (CDM) models. ^^^^^^^ ^ ^j^gp imaging survey in the sky 

More recently, Venemansetal. (2002) have discovered a proto- ^^.^^ j^e Subai'u Deep Field (SDF; centered at 

' Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. 
^ Department of Astronomy, School of Science, University of Tokyo, Tokyo 113-0033, Japan; shimasaku@astron.s.u-tokyo. ac.jp 
^ Research Center for the Early Universe, School of Science, University of Tokyo, Tokyo 11 3-0033, Japan 

National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 
' Institute of Astronomy, School of Science, University of Tokyo, Mitaka, Tokyo 181-0015, Japan 

* Subaru Telescope, National Astronomical Observatory of Japan, 650 N. A'ohoku Place, Hilo, HI 96720, USA 
' Department of Mathematical and Physical Sciences, Faculty of Science, Japan Women's University, Tokyo 1 1 2-868 1 , Japan 

* Research Center for Neutrino Science, Graduate School of Science, Tohoku University, Aramaki, Aoba, Sendai 980-8578, Japan 
' Institute for Cosmic Ray Research, University of Tokyo, Kashiwa, Chiba 277-8582, Japan 

The Graduate University for Advanced Studies (SOKENDAI), Shonan Village, Hayama, Kanagawa 240-0193, Japan 
' ' Department of Physics, University of Durham, South Road, Durham DHl 3LE, UK 
Department of Astronomy, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan 

Astronomical Institute, Graduate School of Science, Tohoku University, Aramaki, Aoba, Sendai 980-8578, Japan 

1 



2 



Shimasaku et al. 



(13''24'"21.'4,+27°29'23") [J2000.0]; Maihara et al. 2001) in 
the R and i bands and a narrow -band filter centered at 7126 A, 
NB71 1, with the prime-focus camera (Suprime-Cam; Miyazaki 
et al. 2002) on Subaru in Mai'ch - June 2001 and May 2002. 
The FWHM of the NB711 filter is 73 A'"*, giving a survey 
depth for LAEs along the sightline of Az = 0.06, or equiva- 
lently 33hjQ Mpc. We observed two field-of-views (FoVs) of 
Suprime-Cam allowing for a wide overlap: the central FoV and 
the northern FoV 
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seeing sizes of the final images are 0."90. 

Object detection and photometry are made using SExtractor 
version 2.1.6 (Bertin & Arnouts 1996). The NB711-band im- 
age is chosen to detect objects. In this Letter, to ensure secure 
photometry and a low contamination for our LAE sample, we 
confine ourselves to objects brighter than NB71 1=25.5 (their 
number is 34,653). We apply the following criteria (Ouchi 
et al. 2003a) to these objects, to isolate LAEs at z = 4.86: 
/;/-NB711 > 0.8, R-i> 0.5, and /'-NB711 > 0, where 
Ri =(R + /') /2. The second criterion reduces contamination by 
foreground galaxies whose emission lines other than Lyman a 
happen to enter the NB711 band (see Fig. 1). The number 
of objects passing the above criteria is 43. Figure 1 plots on 
the Ri-NBlll vs R-i' plane all objects with NB711 < 25.5 
(right) and model galaxies and Galactic stars (left). This figure 
demonstrates that the second criterion efficiently removes low-z 
objects but selects LAEs at z = 4.86 (Ouchi et al. 2003a). 
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Fig. 1. — Two-color diagrams for continuum color ) and nan"ow-band excess color 
(/^/-NB711), where Ri = (R + i) /2 is a continuum magnitude, (left panel) Tracks of model 
galaxies at different redshifts. Red lines indicate model LAEs which are a composite spec- 
trum of a 0.03 Gyr single burst model galaxy (GISSEL96; Bruzual & Chariot 1993) and a 
Lyman a emission (£Wrest(LyQ)=22 A); from the left to right, three different amplitudes 
of IGM absorption are applied: O.Srefr, Tefr, and l.Srcrr, where T^[f is Madau's (1995) me- 
dian opacity. The narrow-band excess in each of the peaks in the red lines indicates the 
Lyman a emission of LAEs 31 z = 4.86. Green lines show six templates of nearby starburst 
galaxies, with different dust extinction {E(B — V) = 0.0—0.7), given in Kinney et al. (1996) 
redshifted up to c = 1.2. The narrow-band excess peaks in the green lines correspond to the 
emission hne of Ha (z = 0.08), [Oin](z = 0.4), H/3(z = 0.5), or [OII](z = 0.9). Two blue 
lines indicate the tracks for model LBGs without Lyman a emission, with O.Srefr and Tefr. 
respectively. Black lines show colors of typical elliptical, spiral, and irregular galaxies of 
Coleman, Wu, & Weedman (1980) which are redshifted from z = to z = 2. Yellow star 
marks show 175 Galactic stars given by Gunn & Stryker (1983), and 24 Kurucz model 
stars, (right panel) Colors of the detected objects. Large circles indicate 43 LAE candi- 
dates, among which red ones have spectroscopic observations, while small dots are for the 
other objects. An open circle with error bars indicates the median colors and their errors 
for the 43 LAE candidates. Our LAE selection criteria are outlined by a pink line. 

The data of the central FoV have already been reduced 
and used to study general properties of LAEs by Ouchi et al. 
(2003a). For the northern FoV, individual CCD data were re- 
duced and combined using IRAF and our own data reduction 
software (Yagi et al. 2002). For the overlapped region, we use 
the data of the central FoV. The total exposure time and the lim- 
iting magnitude (3a on a l."8 aperture) are: 90 min and 27.1 
mag (R), 138 min and 26.9 mag (/'), and 162 min and 26.0 mag 
(NB711) for the central FoV; 120 min and 27.5 mag (R), 110 
min and 27.3 mag (/'), and 144 min and 26.1 mag (NB711) 
for the northern FoV. All magnitudes are AB magnitudes. The 

The central wavelength and the FWHM of the NB711 filter are the values measured for the F/1.86 prime focus. 

" We first observed the central FoV and found in a northern part a large-scale overdense region of LAE candidates. We then imaged the northern FoV with an overlap 
of 15' to check the rehabihty of this overdense region as well as to trace a northern extension of the overdense region. We found consistently an overdense region in 
the data of the northern FoV. The match of LAE candidates in the overlapped region between the two FoV data is ~ 60%, which can be regarded as an estimate of the 
completeness of our LAE detection. 



Fig. 2. — Spectra of five candidates over the range 6900-7300A. Flux of objects A, B, C, 
and D, has been offset for clarity by 4 X 10"'*, 3 X 10"'*, 2 X 10"'*, and 1 X 10"'* ergs 
s"' cm"^ A"' , respectively. The (relative) spectrum of the night sky emission is shown by 
a dotted line; emission lines are not present around 7100A. 



2.2. Spectroscopy 

We selected a b'4> field of the highest overdensity (see §3.1) 
for multi-slit spectroscopy, and made spectroscopic observa- 
tions for five candidates in this field with FOCAS (Kashikawa 
et al. 2002) on Subaru on June 6, 2002, to check our pho- 
tometric selection. Table 1 summarizes their properties. The 
exposure time was 2 hr for each object. The seeing size was 
0."7-0."8. We used the 300 line/mm grating with a dispersion 
of 1 .4 A pixel"' and a wavelength coverage of 4700-9400 A. 
We adopted a slit width of 0."8, which gave a spectral resolu- 
tion of 9.8 A. The continuum flux limit of our spectroscopy was 
6.3 X 10"''' erg s"' cm'^ A"' {5a). 

Figure 2 shows spectra of the five candidates after two-pixel 
smoothing. All the spectra have an emission line near 7100A. 
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Possible candidates (other than Lyman a) for this emission line 
are Ha, H/3, [OIII]A5007, and [OnjA3727. We find that the 
lines are not [OIII]A5007 atz^ 0.4 because of the lack of the 
H/3 line at the corresponding wavelength; from similar reasons, 
they are not either Ha or H/3. Thus, the five objects are either 
an LAE or an [Oil] emitter. 

Table 1 . Five candidates with spectroscopic observations. 



ID 




i'ia) 


NB7II*"^ 


z 




fiO 




(mag) 


(mag) 


(mag) 




(A) 




A 


27.18 


26.32 


24.24 


4.850 


51 


3.5 


B 


28.25 


26.82 


25.33 


4.834 


15 


0.77 


C 


27.22 


26.56 


24.88 


4.856 


14 


0.93 


D 


27.21 


25.96 


25.33 


4.865 


6 


0.49 


E 


27.62 


27.05 


25.26 


4.825<* 


18 


I.I 



(a) AB magnitudes. 

(b) Rest-frame eqttivalent width (A) estimated from the spectrum, assuming the Une to 
he Lyman a. Due to very weak continuum emissions and a relatively short exposure of 
spectroscopy, EW values should have large errors. 

(c) Line flux in units of 10"'^ erg s"' cm"^. We have not examined the profiles of Lyman 
a emissions because of low S/N ratios of the spectra, although A and possibly B appear to 
have wide velocity widths, which may suggest starburst winds. 

(d) z = 0.90 if the line is [On], 



and B shows an asymmetric shape with a blueward cutoff; this 
is a common feature of high-redshift Lyman a emission (e.g., 
Dey et al. 1998). Second, although we cannot rule out the 
possibility of D being an [Oil] emitter on the basis of its spec- 
trum alone, we have seen in Figure 2 that its R - i' color is red 
enough to match the expectation for LAEs. Moreover, Ouchi et 
al. (2003b) find that C and D satisfy the selection criteria for 
LBGs at z = 4.7 ± 0.5 in the V - i' vs i'-z' plane Indeed, 
the V-i' colors oi A,B,C,D are ~ 1.7-2.5, being consistent 
with the color expected for z ^ 5 galaxies; on the other hand, 
the color of £ is ~ 0.5. From the above discussion, we estimate 
the contamination of our sample to be ~ 1/5 = 20% 

As an independent check, we have examined the distribution 
of objects with Ri -NBlll > 0.8 but R-i < 0.4, which are 
probably low-z interlopers, and have found that their distribu- 
tion on the sky is almost uniform and does not correlate with 
that of the LAE candidates. 
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Fig. 4. — Number of 5 > 2 spheres in our survey volume predicted by three CDM mod- 
els. Solid, dotted, and dashed lines indicate the A-dominated, open, and standard CDM 
models, respectively. 

3. RESULTS AND DISCUSSION 

3.1. Large-Scale Structure of LAEs 

Figure 3 shows the sky distribution of photometrically se- 
lected 43 LAE candidates. We find a remarkable, large-scale 
clustering; an overdense region Ues in the ESE - WNW di- 
rection and there are few objects outside this region. To 
quantify overdensity, we estimate the local surface density of 
LAEs, S(x,>')^and compute the surface overdensity Ssix,y) = 
[S(x,)')-S]/_E, where S is the mean surface density of LAEs. 
We adopt as E the mean surface density in our image, although 
a larger area is desirable to obtain a more accurate (i.e., global) 
estimate of E. 

The overdensity contours are drawn in Fig. 3. The projected 
size of the overdense region (S-e > 0) is found to be about 20 
/z^qMpc in width and larger than 50 hj^Mpc in length (comov- 
ing units) since the region seems to continue outside the im- 
age at either side. This elongated overdense structure may be 
a cross section of a wall-like structure which extends along the 
sightline (Note that the survey depth is 33hjQ Mpc). Another 
possibility may be that we are seeing a 'redder' part of a struc- 
ture which is centered at a redshift smaller than that correspond- 
ing to the center of NB71 l(z = 4.86), since all five objects with 
spectroscopy have z ^ 4.86. This elongated structure includes 
a region of > 2 which is approximated well by a circle of 

The five candidates are located in the central part of the SDF, where we obtained deep B,V,z' images as well in March - June 2001 (Ouchi et al. 2003b). 

This value is lower than that estimated by Ouchi et al. (2003a) for a deeper sample, ~ 40%, on the basis of photometric properties alone. This is probably because 
the photometric errors in our sample are smaller than those in Ouchi et al. 
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Fig. 3. — Sky distribution of 43 LAE candidates. North is up and east is to the left. Areas 
of relatively poor quality have been trimmed. A bright star in the upper left comer, around 
which detection of LAEs is impossible, has been masked. LAE candidates are shown by 
circles. Brighter candidates are shown by larger symbols. Dotted, dashed, and solid lines 
correspond to contours of (Se = 0, 1 , and 2, respectively (A top-hat smoothing of Sh^Q Mpc 
radius is made over the LAE distribution to compute the local overdensity). The region 
with 6s > 2 can be approxunated by a circle of i2h^l Mpc radius. 

We conclude from the following reasons that three (A, B, D) 
are convincing LAEs and C is a probable LAE while E may 
be an [Oil] emitter at z ~ 0.9. First, the emission line of A 
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12 /lygMpc radius. Since the minimum mass of this circular re- 
gion is computed to be ~ (47r/3)/9ol2^ ~ 3 x lO^^^yg^© (po is 
the mean matter density of the universe), it may become later 
a massive cluster of galaxies after collapsing to a size of a few 
Mpc. In the present-day universe, clusters are often embed- 
ded in large-scale structures. From these features, it is very 
likely that this elongated region is a proto large-scale structure 
atz~ 5. 

3.2. Implications on Cold Dark Matter Models 

We examine in a simple manner whether or not CDM models 
can reproduce the observed structure. We focus on the circular 
region of > 2, since properties of circular regions are pre- 
dicted easily. We make a reasonable assumption that this region 
is a sphere of 12 /z^QMpc radius in 3-dimensional space with a 
spatial overdensity of 5 = 2. We consider three CDM models: 
a A-dominated model with fio = 0.3, Aq = 0.7, Hq = 70 km s"^ 
Mpc"', (78 = 0.9; an open model with fio = 0.3, Aq = 0, Hq = 70 
km s~' Mpc~\ (Tg = 0.9; the standard model with = 1. Aq = 0, 
Ho = 50 km s~' Mpc"', a% = 0.5. We adopt these values from 
the (Tg - r^o relation derived by Eke, Cole, & Frenk (1996) for 
local X-ray clusters. 

Using Unear perturbation theory, we calculate the number of 
spheres within our observed volume above ^ = 2 as (cf . Steidel 
et al. 1998): 



1 



5l 



Ina, 



exp 



2a? 



dx, 



(1) 



where Vq is the volume surveyed by our observation (1.4 x 
lO^/z^o Mpc^), Vs is the volume of an unperturbed sphere whose 
mass is equal to (l+2/b)^p()(12h'^li)^, (5l is the linearly ex- 
trapolated mass overdensity estimated from the observed mass 
overdensity (we use a fitting-formula in Bemardeau 1994), 
is the rms fluctuation of mass overdensity at z = 4.86 on top-hat 
filter on a scale of the unperturbed sphere, and b is the linear 
bias parameter for LAEs. 

Figure 4 plots N against b for the three models. We find that 
the best- fit value of b that produces the observed number (=1) 
of spheres is ~ 6 (A-dominated model), ~ 4 (open), and ~ 15 
(Oo = !)• We can also estimate the lower and upper limits of 
b by imposing 0.05 < AT < 4.7 (95% C.L. assuming Poisson 



statistics); we obtain b ^ 3 - 16, b ^ 2 - 12, and b ^ 8-37 for 
the A-dominated, open, and fio = 1 models. Thus, LAEs atz — 5 
must be highly biased tracers of mass. If we use the linear bias 
model of Mo & White (1996) which relates the mass of a dark 
halo with its b, we obtain, from the best-fit b values, the mass 
of dark haloes hosting our LAEs to be ~ 1 x 10'^, 2 x 10'^, and 
2 X IO'^/jyqMq for the A-dominated, open, and Oq = 1 models, 
respectively. 

The best-fit values of b we find for our LAE candidates are 
similar to or possibly higher than those for LBGs at j ~ 3-4 
obtained in the following work. From an analysis of a galaxy 
overdensity of ~ lO'^M© at z = 3.09, Steidel et al. (1998) have 
estimated b ^ 4 for the JIq = 0.3 flat model. Analyses of the 
two-point correlation function have given ~ 2 - 3 for LBGs at 
2-^3 and 4 for an flo = 0.3 flat model with the same normal- 
ization as we adopt (Porciani & GiavaUsco 2002; Ouchi et al. 
2001). Higher at z « 5 would be reasonable if we recall that 
b for a given halo mass is predicted to increase with z- Indeed, 
Ouchi et al. (2003b) obtain ^ ~ 6 for LBGs at z ^ 5, compara- 
ble to our value. 

Steidel et al. (1998) and Porciani & GiavaUsco (2002) have 
inferred the mass of haloes hosting LBGs to be of the order of 
IO'^Mq. The b value for z ^ 5 LBGs estimated by Ouchi et al. 
(2003b) also gives a similar mass. Hence, haloes hosting z — 5 
LAEs are roughly as massive as those of LBGs at z ~ 3 - 5 . Pas- 
carelle et al. (1996b) have argued that LAEs (found at z ^ 2.4) 
are sub-galactic objects from their small sizes. Those LAEs 
have also been suggested to be very young (Keel et al. 2002). If 
LAEs at z — 5 are also sub-galactic and if LBGs have relatively 
large stellar masses as suggested by, e.g., Papovich, Dickinson, 
& Ferguson (2001) and Shapley et al. (2001), then the high 
mass of haloes hosting z — 5 LAEs found here will imply that 
dark haloes of a given mass host galaxies of a wide range of 
stellar mass or a wide range of evolutionary stages. 

We would like to thank the referee, William Keel, for use- 
ful comments. We are grateful to the Subaru Telescope staff 
for their invaluable help in our observations. H. Furusawa, F. 
Nakata, and M. Ouchi acknowledge support from the Japan So- 
ciety for the Promotion of Science (JSPS) through JSPS Re- 
search Fellowships for Young Scientists. 
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